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Abstract. The local space density of galaxies as a function of their basic 
structural parameters -like luminosity, surface brightness and scalesize- 
is still poorly known. Our poor knowledge is mainly the result of strong 
selection biases against low surface brightness and small scalesize galax- 
ies in any optically selected sample. We show that in order to correct 
for selection biases one has to obtain accurate surface photometry and 
distance estimates for a large (^1000) sample of galaxies. We derive bi- 
variate space density distributions in the (scalesize, surface brightness)- 
plane and the (luminosity, scalesize)-plane for a sample of ~1000 local 
Sb-Sdm spiral galaxies. We present a parameterization of these bivari- 
ate distributions, based on a Schechter type luminosity function and a 
log-normal scalesize distribution at a given luminosity. We show how 
surface brightness limits and (l+z)^ cosmological redshift dimming can 
influence interpretation of luminosity function determinations and deep 
galaxy counts. 



1. Introduction 

Knowing the space density of galaxies as function of their structural parameters 
(luminosity, surface brightness (SB) and scalesize) is important when: 

1) making comparisons between different galaxy samples, because selection 
functions of extended resolved objects depend on at least two structural param- 
eters. This becomes particularly relevant when comparing samples at different 
redshifts, where (l+z)^ redshift dimming can give rise to strong SB biases. 

2) testing galaxy formation and evolution models, as any successful galaxy 
formation theory will have to be able to explain the spread in structural param- 
eters and their relative frequency in the local galaxy population. 

Many papers have been devoted to the determination of the space density 
of galaxies as function of their luminosity, i.e. the galaxy luminosity function (for 
a recent review see Ellis 1997). In many of these papers one has conveniently 
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ignored the possibility of strong SB biases. Determinations of scalesize distri- 
butions have been scarce (some notable exceptions van der Kruit 1987, Hudson 
k. Lynden-Bell 1991, Sodre & Lahav 1993, de Jong 1995) and often diameter 
distributions are calculated. Diameters distributions are not very useful in sam- 
ple comparisons, as diameters have to be measured at a certain SB level, which 
might differ from sample to sample. Realize for instance, that there may be many 
galaxies that do not have a -025, because there SB is below 25 i?-mag arcscc^^. 

Since the classical paper of Freeman (1970), many papers have been devoted 
to the distribution of SB of disks in spiral galaxies (for review see Impey & 
Bothun 1997). Freeman found that 28 galaxies in his incomplete sample of 36 
had disk central SB values of 21.65ib0.3 i?-magarcscc~^ (sec his review in these 
proceedings). Disney (1976) showed that the limited range in disk central SB 
values might be the result of selection biases. Since then several authors have 
argued that there seems to be indeed an upper limit in the SB distribution near 
Freeman's value, but that the distribution stays nearly flat when going to lower 
SB (e.g. McGaugh et al. 1995; de Jong 1995). Recently this picture has been 
challenged by Tully & Vcrhcijcn (1997), who argued that the SB distribution is 
bimodal, based on K-band data of ~60 galaxies in the Ursa Major cluster. 

In this paper we show that one should not try to separate the distributions 
of luminosity, SB and scalesize, but combine two of these to make bivariatc 
distributions, as any sample will have selection biases in at least two structural 
parameters. 

2. Correcting for selection bias 

Many methods have been devised to correct observed frequencies of object prop- 
erties for selection bias in order to obtain true space density distributions. We 
will here concentrate on the Knax method, where each object gets a weight pro- 
portional to the inverse of its maximum sample inclusion volume (Felten 1976). 
This metod is only correct if the objects are distributed homogeneously in space, 
and therefore the smallest objects in the sample should be visible at distances 
greater than the largest large scale structures in the universe. Homogeneity and 
completeness can be checked with the V/Vmax method (e.g. van der Kruit 1987 
and references therein). Accurate T^nax values can be derived for each object for 
the modern surveys with automatic detection algorithms on digitized data. Each 
object should be artificially blue- or redshifted and be Monte Carlo replaced at 
many positions in the original data set. The recovery fraction of the automatic 
detection routine supplies the volume searched at each redshift shell and pro- 
vides information for confusion limits and Malmquist bias at the survey limits. 
For samples not selected by an automated routine from digitized data (e.g. eye- 
selected from photographic plates), we just have to assume that the selection 
criteria are well behaved when we imagine moving a galaxy in distance. 

Moving more specifically to the distribution of structural parameters of spi- 
ral galaxies, we will use the case of perfect exponential disks in a diameter limited 
sample. More generalised descriptions can be found in Disney & Phillipps (1983) 
and McGaugh et al. (1995). For an exponential disk with physical scalelength 
h and central SB jiQ we find for the maximum distance at which a galaxy can 



2 



lie before dropping out of the sample 

C^max OC (/iUm — /io) /l/6'limi (1) 

with ^lini the sample angular diameter limit measured at SB limit /xiim. As the 
volume where a galaxy is visible goes as d^axi this shows the strong selection 
bias against small scalesize and low SB galaxies. The scalesizes of spiral galaxies 
vary easily by a factor of 10 (de Jong 1996). Therefore, if all scalesizes were 
equally abundant at a given SB, we would have a 1000 times more of the largest 
scalesize galaxies than the smallest scalesize galaxies in a diameter limited sam- 
ple. Luckily nature has not been that cruel to us and there are many more small 
galaxies then large ones. In the case of the SB distribution we have not been so 
lucky, as the SB distribution stays rather constant -at a given scalelength- go- 
ing to lower SB values. Equation (||) shows that, at fixed scalelength, the visible 
volume of a galaxy 1 mag above the SB limit is 125 times smaller than that of a 
galaxy 5 mag above the SB limit. In order to have some number statistics close 
to the selection limit, we had better observe hundreds of galaxies to determine 
a SB distribution. Because we do not a priori know whether the scalesize and 
SB distributions are uncorrelated, we had better make sure that we determine 
the SB distribution at different scalesizes, and so we need at least 1000 galaxies. 

SB measurements are distance independent (at least on local scales); a 
property that sometimes has been used to argue that one can determine SB dis- 
tributions without knowing distances. If the distribution of h is the same at each 
SB level, the /i/^iim factor in Eq. (|l|) cancels out on average and one can make 
relative volume corrections without having to know physical scalesizes/distances. 
Likewise, using total magnitude of an exponential disk M oc /iQ — 5 log(/i), Eq. (|^) 
can be rewritten as 

dmax « (Wim - /io) 10-° '(*'-^«^ (2) 

Again, assuming the SB distribution is the same for each luminosity, one can 
make relative volume corrections (very different from Eq. (||)!) to calculate a SB 
distribution without knowing distances. There is no reason for the SB distribu- 
tion to be independent of either scalesize or luminosity (and we will show this 
is indeed not the case) and therefore Eq. (|l|) & (|2|) show that we need to know 
the distribution of at least one other distance dependent structural parameter 
to determine the SB distribution of galaxies. The reverse is also true: to mea- 
sure the distribution of scalesizes or luminosities we also need to determine the 
distribution of one of the other structural parameters. In order to do so we will 
need surface photometry and distances for a sample of at least ~1000 galaxies. 

In this paper we will use the effective radius (rg, the radius enclosing half 
of the total light of the galaxy) and the average effective SB within this radius 
(</i>e) instead of the more conventional parameters for disks, scalelength and 
central SB. Using the effective parameters has the virtue that one does not have 
to make assumptions about the light distribution in the galaxy (all galaxies 
have an r^, even irregular ones) and avoids complicated bulge/disk decomposi- 
tion issues. The distributions presented here have also been calculated for disk 
parameters alone with very similar results, because most of the objects are of 
late spiral type with insignificant bulge contributions. 
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3. Local space density distributions 

As described in the previous section, one needs accurate surface piiotometry 
and distance estimates for a sample of at least 1000 galaxies to create bivariate 
distributions. The galaxies should be selected in a well defined, reproducible 
and complete way. Data sets obeying all of these criteria are not available at the 
moment, but fortunately peculiar motion studies have produced large data sets 
with accurate photometry and redshifts. We have used the Mathewson, Ford 
& Buchorn (1992, 1996, MFB hereafter) data set, which was selected from the 
ESO-Uppsala catalog, a catalog with galaxies selected and classified by eye from 
photographic plates. We reselected a sample close to the MFB criteria from 
the ESO-Uppsala catalog, allowing us to evaluate incompleteness in the MFB 
sample (some galaxies were not observed due to foreground stars or inability 
to obtain a velocity width). We selected all galaxies from the ESO-Uppsala 
catalog with type 3<T<8 (Sb-Sdm), angular diameter 1.7'<0maj^5', axis ratio 
0.174<b/a<0.776 and galactic latitude |5|>11°. This resulted in a sample of 
1007 galaxies, of which about 850 have /-band surface photometry and redshifts. 

The luminosity, and <n>e values of the galaxies were derived from the 
luminosity profiles and corrected for Galactic foreground extinction using the 
prescription of Schlegel et al. (1998). Corrections for inclination and internal ex- 
tinction were performed following a method similar to Byun (1992). Distance es- 
timates were obtained from the Mark III catalog (Willick et al. 1997) if available, 
otherwise computed from the Hubble distance, with Hq = 65kms~^ Mpc^^. 

Using the Fmax method described in the previous section, we have calculated 
the bivariate density distribution in the (re,</i>e)-plane, which is presented on 
a logarithmic scale in Fig. |^. The paucity of galaxies in the top-right corner of 
the diagram is real, large, high SB galaxies are readily visible. To the bottom- 
left of the indicated 20 Mpc visibility line we are hit by low number statistics; for 
such small, low SB galaxies we are sampling too small a volume to have reliable 
statistics. The distribution shows a dramatic increase in galaxy space density 
going to smaller scalesizes. At a given scalesize, the SB shows a broad distri- 
bution, peaking at about </i>e=21.5/-magarcsec~^. There is some indication 
that the peak in the distribution shifts to lower SB at smaller scalesizes. 

4. Parametrization of the distributions 

In this section we will define a parametrization of the bivariate distributions, 
as an aid to compare distributions derived from differently selected samples or 
to study redshift evolution. We will follow the most simple form of the Fall &: 
Efstathiou (1980) disk galaxy formation theory to derive such a parametrization 
(for extended versions of the theory see e.g. van der Kruit 1987; Dalcanton et 
al. 1997; Mo et al. 1998; van den Bosch 1998). Galaxies form in this theory in 
hierarchically merging Dark Matter (DM) halos, giving rise to a distribution of 
DM halo masses described by the Press & Schechter (1974) theory, which formed 
the inspiration for the Schechter (1976) luminosity function (LF). We will use a 
Schechter LF to describe the luminosity dimension of our distribution function. 

In the Fall & Efstathiou (1980) model, the scalesize of a galaxy is determined 
by its angular momentum, which is acquired by tidal toques from neighbouring 
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Figure 1. The space density distribution of Sb-Sdm galaxies as 
function of effective radius and average SB within that radius. Galaxies 
with exponential disk, having structural properties indicated by the 
line, can be seen out to 20 Mpc before dropping out of the sample. 

DM halos in the expanding universe. The total angular momentum of the system 
is usually expressed in terms of the dimensionaless spin parameter (Peebles 1969) 

\ = j\E\^iHi;^'^G-\ (3) 

with J the total angular momentum, E the total energy and Mtot the total mass 
of the system, all of which are dominated by the DM halo. N-body simulations 
(e.g. Warren et al. 1992) show that the distribution of A values acquired from 
tidal torques in an expanding universe can be well be approximated by a log- 
normal distribution with a dispersion a\ ~ 0.5 in ln(A). 

A few simplifying approximations allow us to relate each of the factors in 
Eq. (^) to our observed bivariate distribution parameters. A perfect exponential 
disk of effective size rg, mass M^, rotating with a flat rotation curve of velocity 
Vc has Jd oc MdreVc- We assume that the specific angular momentum of the 
disk is equal to the specific angular momentum of the dark halo. From the virial 
theorem we get E oc V^M^ot- If we assume that light traces disk mass {M^ oc L) 
and that disk mass is proportional to total mass (Mtot oc -M^d), we only need the 
Tully & Fisher (1977) relation (L oc V/, with [3 ~ 3 in the I-passband) to link 
the spin parameter A to our observed bivariate distribution parameters. 

These approximations give A oc reL~^'^/^~^^ ~ rg-L"^/^. As A is expected to 
have a log-normal behavior, this means that, at a given luminosity, we expect the 
distribution of scalesizes to be log-normal, and that the peak in the r^ distribution 
shifts with ~ L~^/'^. This is exactly the behavior that is shown in Fig. ||, where 
the function over-plotted on the data shows the log-normal behavior at each 
luminosity bin, shifting by L~^^^ between the luminosity bins and where the 
height is determined by the Schechter LF. 
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Figure 2. The space density distribution of effective scalesizes in 
different bins of total /-band luminosity as marked in the top-right of 
each panel. The histograms represent the data with errorbars showing 
the 95% confidence limits due to distance and Poisson errors. The 95% 
confidence upper limits were calculated using exponential disks and the 
survey limits determined from the photometry. The dashed line shows 
the bivariate distribution function described in the text. 



The function plotted is the result of the well known x-by-eye fitting method, 
and the detailed parameters will definitely change when a full fitting technique 
has been developed that takes the Poisson errors on the data points into account. 
For reference we list here the full bivariate function in magnitudes, and the 
parameter values giving a good approximation to the data: 

cf^f '^0 f l^ log re/re. - 0.4(M - M,)(2//3 - 

^{re,M)dlogredM = -=exp -- 

crA\/27r 2 o-A/ln(10) 

lQ-o.4{M-M,)(a+i) exp(-10-0-<^^-^^*)) dlogr.dM, 

with the first line representing the log-normal scalesize distribution and the 
second line the Schechter LF in magnitudes (M). The x-by-eye parameters are: 

$0 = 0.002 Mpc"^ a = -1.25 /3 = 3.0 (slope Tully-Fisher relation) 

M^, = —22.3 /-mag r^* = 6.7 kpc ax = 0.3 

The width of the spin parameter distribution {a\) we need is less than what is 
typically found in N-body simulations. 
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Figure 3. Luminosity functions computed by integrating the bivari- 
ate distribution function down to the indicated I-band central surface 
brightness hmits. The uppermost hue shows the integrated total LF. 

5. Discussion & conclusions 

Using our parameterization we can now estimate the effects of SB limits on 
determinations of LFs and on deep galaxy counts, especially in the context of 
(l+z)^ redshift dimming. If we were to determine the LF from a galaxy sample 
with low SB galaxies cut out, we would underestimate the faint end of the LF, as 
most low SB galaxies are also low luminosity systems. SB cuts are relevant for 
redshift surveys selected from shallow (photographic plate) material (resulting 
in implicit cuts) and for most fiber based redshift surveys which often have 
explicit SB cuts (e.g. the Las Campanas and Sloan Surveys). How SB cuts 
can effect LF determinations is shown in Fig.^, where we have integrated our 
parameterization down to the indicated /-band central SB. 

Figure^ suggests that the SB cuts typically present in local surveys do not 
dramatically effect LF determinations (using for a typical spiral i?-/~1.7), espe- 
cially taking into account that most spirals have some central light enhancement 
due to the bulge, making detections easier. The situation changes however when 
we move to higher redshifts and have to take (1+z)^ cosmological redshift dim- 
ming into account. At z=l our SB limit has already shifted 3 magnitudes up, 
and 6 magnitudes by the time we reach z=3. This means that even for the deep- 
est image available at the moment -the Hubble Deep Field- the SB cut at z=3 
(the [/-band dropouts) runs at about 21 /-mag arcsec^^ (using a K-correction of 
an unevolved Sb galaxy). This limit makes a considerable fraction of galaxies in 
Fig. H undetectable, if we put this local galaxy population unevolved at z=3. 

Tully & Verheijen (1997, these proceedings) have argued that the central SB 
of galaxies shows a bimodal distribution, in particular when looking at /C-band 
data. We do not see such bimodality, independent whether we use their proposed 
bimodal dust extinction correction, we use only the 200 most face-on galaxies 
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with the smaUest extinction correction, we use bulge/disk decomposed param- 
eters or effective parameters. In the many ways we have looked at the MFB 
data set, we have never seen any bimodahty in the SB distributions. Whether 
the bimodal effect is the resuh of the special Ursa Major cluster environent that 
was studied or an unlucky case of low number statistics remains to be seen. 

The simple parametrization presented in this paper gives an accurate repre- 
sentation of the observed bivariate distributions, independently of whether one 
believes in hierarchical galaxy formation models or in CDM-like universes. A 
detailed analysis of galaxy formation in CDM-like universes paying attention to 
bivariate space density distributions will appear in Lacey et al. (1999). 
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